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Abstract 



The Witten Laplacian in one dimension is studied further by methods of resurgent analysis 
in order to approach Fukaya's conjectures relating WKB asymptotics and disc instantons. We 
carry out explicit computations of exponential asymptotic expansions of exponentially small (i.e. 
0(e^^/'*), c > 0, /i — 0+) eigenvalues and of corresponding eigenfunctions of the Witten Laplacian; 
a general algorithm as well as two examples are discussed. 

1 Introduction 

We are continuing the project started in |G11] . where we proposed to study the Witten Laplacian 
by methods of resurgent analysis in order to prove conjectures by Fukaya ,1-0."^ §5.2] relating WKB 
asymptotics and disc instantons. The reader is referred to the introductory section of [Gllj for 
philosophy and motivation, as well as for a brief review of resurgent analysis. 

The present paper is purely computational; it relies on techniques of resurgent analysis and the 
theory of exponential asymptotics whose rigorous mathematical justification is a subject of current 
research. 

Resurgent functions and exponential asymptotics. We will continue to use the termi- 
nology of [Gllj which is a blend of |CNPj and jShStj . In particular, by a resurgent function of the 
semiclassical parameter h — )• 0+ we will mean an equivalence class of holomorphic functions (f){h) 
representable as a Laplace transform 



of an endlessly continuable function <1*(^) ("major of (p") along an infinite contour 7. See |G1H 
Sec. 2] for the precise definition of the equivalence relation used, the meaning of the words "endlessly 
continuable" , and the choice of 7. 



m 
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Under favorable conditions on singularities of ^*(^), a resurgent function </>(/i) can be developed 
into an exponential asymptotic series 

^ Cj 

4>{h) ~ e~"h (aoj + + a2j^^ + •••)• (1) 
i 

In the present paper we will also see resurgent functions where the RHS of ([T]) contains series in h 
and /iln h. 

If (/>(/i) depends on any additional parameters, then so does ^{i), and the RHS may exhibit 
apparent discontinuities; this is known as a Stokes phenomenon, cf. |Glll Sec. 3]. 

We think of resurgent analysis as a way of making mathematical sense of exponential asymptotics 
and describing the Stokes phenomenon. In the present paper we will work mostly with expansions as 
m the RHS of and use a description of the Stokes phenomena available from the previous literature; 
we will cite analytic statements without discussing their current rigor status. 

Spectrum of the Witten Laplacian. Consider a generic enough, in the sense of |G11| . real 
trigonometric polynomial f{q) G M[sin27rg, cos27rq] with n real local minima and n real local maxima 
on [0, 1), and associate to it the following /i-differential operator called the Witten Laplacian 

P = -h''dl + {ff-hf", h^0 + . (2) 

Motivated by Morse theory, we are interested in exponentially small, or low-lying eigenvalues E{h) of 
P with periodic boundary conditions, namely those eigenvalues that can be estimated < E{h) < e~s 
for h — )• 0+ and c > 0. 

In |G11| we have shown, modulo standard black boxes in resurgent analysis, that P has n 
exponentially small resurgent eigenvalues E}^[h) and that the corresponding eigenfunctions (pk{Q,h) 
are resurgent with respect to h for q There we also sketched a method for calculating 



Ek{h) which we recall in section 2.1 For technical reasons, in [Gllj we worked with the semiclassical 
parameter h satisfying < arg/i << 1, \h\ — )• 0. 

Because of the Stokes phenomenon, we cannot hope to write asymptotic expansions of eigen- 
functions of P valid for all q G [0,1); for every line segment between two consecutive critical points 
of f{q) a different exponential asymptotic expansion is valid. Those asymptotic expansions were not 
discussed in |Gllj and will be calculated here. 

Statement of results. This purely computational paper has two main results, both of which 
contribute to our computational dexterity when it comes to explicit exponential asymptotics of eigen- 
values and eigenfunctions of the Witten Laplacian corresponding to a trigonometric polynomial 
f{q) with n real local minima qi, ...,q2n-i and n real local maxima q2,--,q2n on the period, where 
> qi < q2 < ■■■ < q2n-l < q2n < 1- 

The first result is a computation of the connection coefficients for the equation P^ = hEj-ip, 
Er G C, that matches exponential asymptotic expansions of its solutions on the intervals {qj-i, qj) and 
{qj,qj+i) (called connection coefficients across the point qj) to order 0{h). Less precise results and a 
somewhat shortcut treatment of [Gil] was stopping us from confidently proceeding to the calculation 
of periodic eigenfunctions of the Witten Laplacian, see below. 

Theorem 1.1 Let qi be a real local minimum and q2 be a real local maximum of f{q). Then: 



a) with respect to the basis of formal WKB solutions (28) the connection coefficient, defined by (29), 
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across the point qi is expressed by the formula ( 36 ) ; 



b) with respect to the basis of formal WKB solutions (37) the connection coefficient, defined by (38), 



across the point q2 is expressed by the formula (42) 



We warn the reader again that this theorem is proven modulo standard black boxes in the 
complex WKB method; here we present a calculation whose complete justification is a topic for future 
research. 



In finite time, the formulas (36), (42) results could have been made more precise both with 



respect to h and Er, see Sec, 5.6 



The intrigue of this calculation is the conjectured analyticity of the reduced connection coef- 

about the 



ficients {c'^y'^iE, h), cf.(31), and cl'"^{E,h), cf.(40), with respect to E. See remark 



5.1 



current mathematical status of this analyticity statement known as the Sato's conjecture, and sec, 5.5 



for a partial argument in the case of the Witten Laplacian. Thus, our calculation is a nice explicit 
example of the Sato's conjecture. 

The second result consists in explicit calculation of the low-lying eigenvalues and the correspond- 
ing eigenfunction of the Witten Laplacian in two specific examples. 

A bit of notation: let mean terms of exponential asymptotics whose exponential prefactors 
are < es ; let us write (j){h)'^eh for a G M if the leading exponential term in the exponential asymptotic 
expansion of is times an expansion in terms of h and hlnh. 

Example 1 (section^^. If f{q) = [sin27rg + cos47rg], then the operator ([2]) has two exponen- 
tially small eigenvalues 

^(^^=0 and E^'^^ ^ e'^{<o^h + o{h^)) + e-^{ hln h + 0{h^)) + £~ 



25 
87r . 



TT 



The eigenfunction corresponding to E^^'^ = is e""^ , and the eigenfunction corresponding to -E^^) has 
the following asymptotics (actually, more precise information on the asymptotics of the eigenfunction 
is easy to obtain from our calculations in section Isl): 



^ ^ (0' 3 ~ ^arcsini) 



1 • 1 l^ 

- 2^arcsm^, 2, 
,2,| + ^arcsin|; 
'I + ^ arcsin 1, 1) 



4' 



(1 + £~s^)i>+{q, h) + ie^{^ + 0{h))i;_{q, h); 
(1 + £-i^)4>+{q, h) + ie-^{^ + 0{h))^.{q, h); 
(_1 + £-^)^^{q, h) + ie~^{^ + 0{h))i^.{q, h); 
(_1 + h) + ies^ij^ + 0{h))^.{q, h), 



(3) 
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(please keep in mind that 



where the basis ip±{x, h) of formal WKB solutions is introduced on page 

/here(9+ |) = /p|3T](g))- 

Post factum it turned out that we did not need 0(/i)-level contributions to the connection 
coefficients to obtain these formulas; perhaps, this is the case for all sufficiently generic examples of 



f{q). Using the full precision of the formulas (36) and (42) would have given us one more term in the 
(/i. In /i)-expansions multiplying e~8^ and e~8^ in the expression for E^^\ 

Example 2 (section^. If f{q) is a real trigonometric polynomial with two local minima qi,qz 
and two local maxima q'2 , , < (71 < ^2 < 93 < 94 < 1 j such that 



/(9i) = 0, f{q2) 



1 

2' 



b 
2' 



/(94 



a 
2' 
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where 



< 6 < a < ^ and 2a < 3b, 



then the operator ^ has two exponentiahy smah eigenvalues 



i?W=0 and 



b — a 

~e h , 



The eigenfunction corresponding to E^^^ = is e ^ \ and the eigenfunction corresponding to E^'^^ 



has on {qj,qj^i), := qi + I, j = 1,..,4, the exponential asymptotic expansions D 
where the basis ip±{x,h) of formal WKB solutions is introduced on page 



D^£^ are given by (|68 
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and approximations to 



In remarks 8.2[ 9.1 we put our finger on the specific algebraic reason why methods of complex 
WKB resurgent analysis are essential for such a calculation and why they look more powerful than 
C~ methods of, e.g., |HKNn4j . 

In conclusion, the ability to perform explicit calculations developed in this paper will be needed 
in our future work towards Fukaya's conjecture. Remarks 8.2, 9.1 and computations leading to them 
may be of independent pedagogical interest. 

The structure of the paper follows the outline given in section [2] 



2 Method for calculation of the spectrum of the Witten Laplacian. 
2.1 Eigenvalues 

We will now review a method to calculate exponentially small eigenvalues of the Witten Laplacian 
proposed in |G11| . 

The semiclassical parameter h is assumed to satisfy < arg/i << 1, \h\ — t- 0. 

Let f{q) be a real polynomial in sin27rg and cos27r(7, with n real local minima qi, q2n-i and n 
real local maxima q2, ■■, q2n on the period, where < qi < q2 < ■■■ < q2n-i < Q2n < 1- We will assume 
that /'(O) ^ 0. 



Step 1. Formal solutions of the Witten Laplacian. Let E he a, complex number, |£'| sufficiently 
small. We choose two formal WKB solutions 

V^±(i?, q, h) = e^^ia^iE, q) + afiE, q)h + a^{E, q)h^ + ...) (4) 

of the differential equation 

Pi>{E,q,h) = EtP{E,q,h). (5) 

The words "formal WKB solution" mean that iIj± satisfy ^ in the sense of formal power series in h. 
Dependence on E will be often suppressed in our notation. Actually, different choices of the formal 
solutions will be convenient in different chapters, so the notation ^-t (E, q, h) will not be kept beyond 
this SeciSl 
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The equation ([5]) is a stationary Schrodinger equation with the potential [f (q)]"^ — E — hf" (q) . In 
the standard terminology of the WKB analysis, points q for which the /i-independent part [f'{q)]'^ — E 
of the potential vanishes are called turning points of the equation ([s]) . They are called simple or double 
turning points according to the multiplicity of the zero of the function [/'((?)]^ — E. Higher-order 
turning points are also studied in resurgent analysis but they will not be relevant for our discussion. 

When E = 0, the equation ([s]) has 2n double turning points on [0, 1) located at the critical points 
of f{q). When we deform E away from zero and make < \E\ « 1, each of these double turning 
points gives rise to a pair of simple turning points. This is the origin of several delicate features of 
our analysis. In particular, the ingredients S{E,q), a^{E,q) on the RHS of Q are ramified analytic 
functions of E,q. 

Section [3] is devoted to the explicit form of the formal solutions |4] and to their properties as 
multivalued functions of E and q. 



Step 2. Connection problem; limit E = hE^. This next step crucially requires the study of the 
equation ^ for complex values of E and q, but its outcome can be explained in term of q confined to 
the real axis and E >Q. 

The critical points < (/i < ... < q2n < 1 of f{q) are double turning points of the equation ([s]) 
for E = {). For < << 1, they split into pairs of simple turning points < qi {E) < qi{E) < 
q:,{E)<...<q+^{E)<l. 

The content of the Stokes phenomenon is that actual (as opposed to formal) resurgent solutions 
ip{E, q, h) of ([5]) have different asymptotic expansions in terms of the formal WKB solutions il^±{E, q, h) 
on different intervals between the turning points; see |Gllt §3] for definitions. E.g., 

i,{E, q, h) ~ A+{E, h)^+{E, q, h) + A.{E, h)ij_{E, q, h) for g7_^ < g < q.{E), 

and 

i^{E,q,h) - B+{E,h)^P+{E,q,h) + B^iE,h)^^{E,q,h) for q+{E) < q < qJ+.iE); 

here we let q+iE) = q+^{E) - 1 and q^^^^iE) = q^{E) + 1. 

In |Gllj we conceptualize coefficients A±[E)^B±{E) as resurgent symbols, see |G1H §2.2]. For 
< |£^| << 1, they correspond to exponential asymptotics as on the RHS of ^. 

The connection problem consists in finding a 2 x 2 connection matrix C^'^^\E, h) of g-independent 
resurgent symbols such that 

The delicate and, to our knowledge, only partially understood on the rigorous level step is to 
replace a number E by an expression hE^, Er £ C; here we follow the notation of [DDP97| . |DP99j . 
where "r" stands for "reduced". The /i-independent part of the potential (/')^ — hf" — hE^ is now 
(/')^, and so the equation 

h) = hEri^iq, h), Er £ C (6) 

has double turning points at qj, j = l,...,2n. A careful calculation of the limit yields connection 
matrices C^'ii\hEr , h) . 
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We discussed this calculation in |G1H §7], but now we would like to have more accurate results 
- more terms in the asymptotic expansions. This is accomplished in section [5| 



Step 3. Quantization condition. We will now write down conditions on the complex number Er 
for the equation ([g]) to admit a periodic resurgent function solution. This requirement is equivalent 
to the existence of coefficients (resurgent symbols) not simultaneously zero, satisfying 

{^l^+{hEr, q + 1, h),'tP_{hEr, q + 1, h))C^'^''^\hEr, h)....C^'i^\hEr, h) 

= {^P+{hEr,q,h),^l^.{hEr,q,h)) ( ;^+|J;| ) . (7) 

We call the 2x2 matrix 

V 'i/j+{hEr,q,h) 'il^_{hEr,q,h) ) 

the transfer matrix. 

A condition on E^ that there exists a nonzero pair {A^[h), A^{h)) satisfying Q is called the 
quantization condition; it can be written in the form 

det{F{Er,h) - Id) = 0. (9) 

The algebraic structure of the quantization condition was made explicit in [Gilt §8] using the 
specific form of the entries in C^'^i\hEr, h). The quantization condition can be written in term of the 
following 4n + 1 quantities: 

i) ^j{hEr,h), j = l,...,2n - monodromies of one of the formal solutions tp±(hEr,q,h) around the 
turning points qj; 

ii) Tj{hEr,h) for j = l,...,2n - limits of monodromies of one of ij)±{hEr,q,h) along the "tunneling 
cycles" connecting qj and qj+i; 

iii) l + K{Er,h)Er = tp±{hEr,q + I, h)/tp±{hEr,q,h) ioT one oi the choices Tp+ or ip^, c(.(A7). 
We will recall precise statements in sec, 7.1, cf.(48). 



Step 4- Solving quantization conditions by means of a Newton polygon. The equation ([9]) is 
satisfied for Er = 0, but otherwise the requirement that Ej. should be an /i-independent complex 
number is too restrictive. We will now solve the quantization condition ^ allowing Er to be any 
exponentially small resurgent function, Er = Er{h) = 0(e~s), c > 0. 

As in |G1H §9], we write the quantization condition in the form 

T{Er,h):=^E^e^nak4h) = 0, (10) 

k,c 

where k ranges over positive integers and c over a subset of M bounded from above and without 
accumulation points. 
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Figure 1: An example of a Newton polygon representing the quantization condition 



We plot the points (/c, c) corresponding to nonzero afc^c(^) on a coordinate plane, and consider 
the Newton polygon M which is the convex hull of 

[J \k\ +oo) X (-oo; c], 

(fc,c):afc,,(/i)^0 

By the leading term of M we will mean the term E^ehak^dh) corresponding to the leftmost 
vertex of the horizontal edge of AA; such a term necessarily has degree n in Er, \G11\ §10]. In |Glll §9] 
we have written up an iterative procedure for finding exponentially small solutions oi J-{Er{h), h) = 0. 
In a nondegenerate case there will be n such solutions E^^ = 0, Ei'^\ and for k = 2, .., n 

E(k) ^ e-tfc/'^5(/i^ In h) + e.l.o.t., 

where tk is one of the slopes of an edge of J\f, b{h, In h) is an expansion in terms of h and In h, and 
e.l.o.t. means terms of exponentially lower order. 



In sec 8.1 we will perform this calculation in details for a specific example. 



We finally note that substituting E = hEr{h) into the original equation ([9]) and into the resurgent 
solutions ip{E, q, h) gives resurgent eigenfunctions of P corresponding to this resurgent eigenvalues. 
Rigorous results partially justifying this fact are given in |G12j . 

2.2 Eigenfunctions 

Eigenfunctions of the Witten Laplacian corresponding to an eigenvalue hEr{h) can be calculated as 
follows. 

If Ey{h) satisfies ([o]), then the kernel of the matrix F{Er{h), h) — Id, where F is the transfer 
matrix ([s]), is nontrivial. Since this is a 2 x 2 matrix, an element of its kernel can be read from the 
coefficients of the matrix. 
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Clearly, A^(h)'il^+{hEr{h), q, h) + A-{h)ip-{hEr{h), q, h) will be an asymptotic representation of 
the eigenfunction for < g < gi. 

Using the connection formulas to compare asymptotic representations of eigenfunctions on dif- 
ferent intervals {qj,qj+i), we conclude that on an interval {qj,qj+i) the asymptotic representation of 
the same eigenfunction is given by 

{tl;+{hEr,q,h),tl;_{hEr,q,h))C'-'^^\hEr,h)....C'^'''\hEr,h) 

We present results of this kind of computation in two examples, sections |8.2| and [9} 

Remark 2.1 In order to obtain nontrivial results on all intervals {qj,qj^i), we have to work with 
exponential asymptotic expansions of A±{h) beyond the leading exponential order which in turn depend 
on the subleading exponential terms in Er{h). This suggests that an analogous computation would 
be hard to perform by C°° methods, i.e. without complex WKB. 

3 Notation and formal WKB solutions. 

3.1 Notation, cuts, signs, and branches. 

For the purposes of calculation performed in this paper, it is enough to limit our considerations to a 
neighborhood of the real axis in the g-plane. 

Let us recall the notation of jGllj . Let f{q) G ]R[sin(27rg), cos(27r(7)] be a real trigonometric 
polynomial, with n real local minima qi, ...,q2n-i and n real local maxima q2,--,Q2n on the period 
[0, 1), where < gi < ^2 < ••• < q2n-i < q2n < 1- We require /"(qj) / 0. 

In this section [3] we will discuss formal WKB solutions of 

:= [-h^d^g + (/')' - hf"] V' = Ei;, (11) 

where E is a complex number and -0 = i^iE, q, h). 

For E ^ and \E\ sufficiently small, the classical momentum p{q) = ^JE — (/'((?))^ is defined 
on a two sheeted cover of the complex plane of q. For = 0, the two determinations of are 
if'{q), and one can think of the Riemann surface of p{q) as of two separate sheets having contact at 
points qj where f'{qj) = 0. 

Our formulas will be written for E > 0; analytic continuation to other values of E will be 
implicit. 

For every fixed E, the ramification points of p{q) coincide with solutions of the equation E = 
(fiq))"^ which we call turning points of the equation Pip{q,h) = E'ijj{q,h). li E = {), we speak of 
double turning points because {f'{q))'^ has double zeros at the critical points qj, 1 < j < 2n, of f{q). If 
< E << I, the double turning points qj split into pairs of simple turning points qj{E) < qj < q^{E). 

The Riemann surface of p{q) can be described as the plane with cut connecting qJ to q^ and 
going a little below the real axis. To specify the determination of p{q) on the first sheet, we define 
Arg {E — (/')^) for real values of q on figure [2j As — )• 0, on the first sheet ip{q, E) — )• f'{q). 

On our pictures we will draw contours (or their parts) on the first sheet of p{q) as solid curves, 
and contours on the second sheet of p{q) as dashed curves. 
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^1 ArgE-(f)^ = ^2 ArgE-(fr = -271*^2 



Arg E-(f 'f = n \ ^ Arg E-(f 'f = -n 

cut cut 



Figure 2: Choice of Arg E — (/')^ on the first sheet of the Riemann surface of p{q) for a fixed E, 
0<E «1. 

3.2 Formal solutions of Pip = Eip. 

In order to find a formal WKB solution of ( |11[ ), we will be looking for a series 

y{q, h) = yo{q) + hyi{q) + h'^y2{q) + ... (12) 

solving the equation 

(P - E) jexp -J ^E-{nq)Y + y{q' , | = (13) 

in the sense of formal power series. The equivalent condition on y(q, h) is given by the Riccati equation 

2y{q)i^E-{f'Y = -hy{q, hf + i^===, " hy'{q) - f" 

which yields a recursive procedure for calculating yj. In particular, 

ff" f" 



yo{q) 



2{E-im 2i^E-{f'f' 



8z(i; - (/')2)5/2 2(i;-(/')2)2 8i(£; - (/')2)3/2 4i(i; - (/')2)3/2 4(£;-(/')2)' 
etc. 

We have two choices of p{q) = \J E — (/'(g))^ corresponding to the first and the second sheet of 
the Riemann surface of p(g) introduced in Sec]3.1[ Accordingly, we have a two formal solutions 



±(ii;, q, h) = exp |£ '-^E-{f'Y + y{q')d<^^ 



where the sign "+" corresponds to the first sheet and the sign "— " corresponds to the second sheet. 
For definiteness, we put qo = 0, although any other point in the interval {q2n — 1, ^i) would work just 
as well. It is customary to say that (t)±{E, q, h) are normalized in such a way that 4'±{qo) = 1- 

For < E « 1, (j)-^-{E, q, h) becomes (l)-{E,q,h) times a g'-independent factor Cj{E,h), and 
vice versa, when we cross a cut between qj{E) and qj'{E). Further, to make (p±{E,q,h) univalued 
functions of q, let us introduce additional cuts as shown on figjSj 
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+ 
^2 



tt^ 



domain of ^ (E,q) 



+ 

^2 



t/3 
— 



+ 
t?3 



+ 
^4 



domain of ^ (E,q) 



Figure 3: Domains with respect to the q variable of 4'^{E, q, h) and 0_ {E, q, h), < E « 1 
3.3 Formal solutions of Ptp = hEr^. 

Inserting hE,., E^ G C, into the formal solutions (j)±{E,q,h) and regrouping the terms accordingly to 
the powers of h yields formal solutions (t)±{hEj.,q,h) of the equation 

{-h^dl + [/f - /i/")<^ = hEr^. (14) 

Their domains will be taken the same as on figp except cuts between qj and g^" are now shrunk to 
punctures. For the equation Pip = hErtf^, we wm take p{q) = ±f'{q) and speak of two sheets of the 
Riemann surface oi p{q): upper sheet with p{q) = f'{q) and lower sheet with p{q) = —f'{q). 



4 Monodromies of formal solutions. 

For a fixed -E, let p{t), < t < 1, be a path on the Riemann surface of p{q), and let (j){q, h) be a formal 
WKB solution as in Sec )3.2[ but now understood as a multivalued analytic function on the Riemann 
surface oip{q). We call (p{p{l)) / (j){p{0)) the formal monodromy, or the monodromy of a formal solution 
along the path p, and we call Sp{h) satisfying (j){p{l)) / (f){p{0)) = exp[27risp], the monodromy exponent 
along the path p. 

Definition for < << 1 of the paths whose monodromies we would like to compute is given 
on figjl] and figjs} On fig js} e will be taken as a small complex number with < arg e < vr; this 
restriction will be useful later in order to specify which Stokes region the point qk — £ belongs to. 

Notice that 7^ and 7^, also make sense when E = whereas cj^ and cr^ get "pinched" when 
E-^0. 




Figure 4: Paths 7^ and 7^. 



4.1 Some Taylor series 



In this Sec, 4.1 we collect various Taylor series and simple formulas among them that will be relevant 
for the calculations later on. 
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Figure 5: Paths Ufc and o"^. 



Let € M and /'{qi) =0. For q near qg the substitution 

u = -fiq) 

is one-to-one, and thus u can be taken as a local coordinate near qg. 

Let us introduce the numbers aj{qg) (sometimes written as aj if the index d. is clear) by 



-{q - qd = a,{qg)nq) + \ar{qg)[nq)f + \a2{qd[nq)f + 



In particular, 



1 /(^^H?^) f^'\<ld 3[/(3)(g,)] 



2 



'^O — — ^1 «1 — 77777 — 777 5 ^2 



/"(%)' [f"{qdr 2[/"(g,)]4 2[/"(g,)]5 '•- 

In general, the Lagrange inversion formula allows us to write a general expression for a„. 

It follows then that 

1 °° 
- 7//^ = «ofe) + + a2((Z^)[/'(<7)]' + ... = 5]a,(-l)%^ 

i=o 

and 

oo oo 

i=i j=i 
We similarly introduce coefficients bj = bj{qi) by the requirement that 

- /"(?) = hiqi) + + b2iqe)u^ + .... 

should hold near qg. In particular, 

, f/V ^ , /^'Hg^) , [f^'Hqi)?-f"{qe)f^'Kqe) 

We obtain by differentiation 

oo 
j=0 

Finally, for A = —f'{qe — e), we have 

f{qi)-f{qi-e)= / = / nE(-l)%V](in = ^(-1)^- 

J Ql — e JO ■_n r, 



J + 2 

i=o i=o 



^We thank Prof.S.Garoufalidis for this remark. 
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4.2 Formal monodromy along at, k odd. 

In order to have fewer indices in the notation, we wih treat the representative cases of k = 1. 
We have: 

oo 

j=0 

where {AS){E), {Ayj){E) are analytic functions of E at least for small positive E defined by 

ASiE) = z [ - U'Ydq', Ay,{E) = [ y,(q)dq, (19) 

As we will be taking limit for — )■ 0, we will be interested in approximate values of AS{E), 
Ayj{E) for E 0+. 

4.2.1 Calculation of the A5 summand in 27riS(j^. 



We have 



AS= f i^E-{nq))Hq = 2 T [-^{f'Y-E]dq 



(change of variables u = —/'((?), u = \/^ cosht, A = —f'{qi — e); aj = aj{qi) were defined in section 



4.1) 



oo oo „o 

2^(-l)^aj / uWu^-Edu = -2^{-iyaj / E 2+^ coshUsmh^ tdt 

■_n J A ■_n ^arccosh (A/VE) 



j=0 ^ j=0 



arccosh {A/\/E) 



(use (72), (18) and (70) ) 



r n/ \ I, , XT , E E"^ 2 a E , A? E E'^ E'^ ^ , 2A , , 



J = 1 or > 3 



4.2.2 Calculation of the Ayo summand in 27:isa^. 
Recall that 

, , ff" f" 

ym) 



2{E-{rf) 2i^E-{f'Y' 
Lemma 4.1 We have 

i-nq i_-e)) _ TTi 
E 2 ' 



Ayo{E) = / yo{q)dq = arccosh 



where the branch of arccosh is chosen so as to coincide with the principal real value of of arccosh for 
E > and qi — £ on the real axis immediately to the left of q^ . 
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Figure 6: Contour a in the complex plane of u, sec 



4.2.3 



Proof. Integrating the first summand in uq, we have 



1 vri 

= 27ri = . 

aa, 4 2 



i 2(/-W)^^=-^^(^-(^' 
To integrate the second summand, use a substitution u = —f'{q) and A = —f'(q): 

r f'dq f f'dq (-du) A 

J a, 2i^E-{f'Y J a, 2^{f'Y-E J A ^/v?^ Ve' 

(In the second term of this line the arithmetic square root is meant when E > and when q is real 
immediately to the left of q^ .) Subtracting the latter value from the former, we obtain the statement. 
□ 

4.2.3 Calculation of the Ayi summand in 27:isaj- 
Recall that 

. . ^ 5{fTif"r nrr ir? /^^^ 

y^^"" 8i(£; - (/')2)5/2 2(S-(/')2)2 8z(E - (/')2)3/2 Ai{E - {f')^)^/^ A{E - {fY)' ^ ' 

In the integral /^^ yi{Q)dq let us make a substitution u = —f'{q), write —f"{q) = + biu + 
h2V? + ... as in ( [l7| ) , and put A = —f'{qi — e). Choose the contour a in the u-plane as on figjoj if A 
is close to the real axis, in the formulas below arg(i? — v?)^/"^ ~ ^ at the beginning of the path a and 
arg(£' — ti2)i/2 PS — I at the end of the path a. Then 



yi{q)dq 



du. 



Since ^e-u'^]'^ — '~ 4(E-u^]'^^ — ^® ^ ^^^^ differential of a function univalued in u, this part of 
the integrand can be dropped, and so 



yi{qm \\ 8i{E-u^)5/^ 8i{E-u^)V^ At{E-u^)V^^ 
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We prefer to rewrite the denominator in terms of (u^ — E)^^"^ which is positive for u real, u > ^/~E and 
close to the beginning of the path cr, and negative close for u close to the end of a: 



Integrating by parts twice while using formulas (|73[)-([76|), we obtain 



du. 



24(^2 - E)3/2 



bo u 



UE Vn2 - E 



da 



^ (2 - j)bj bju^-^ 

^ 24 (^2 - S)i/2 
j=l > 



da 



2A{v?-EYl^ 

12(^2-^)3/2 ' 6E(A2-£;)l/2 ' 12(^2-^)1/2 " 

tl 
1 



The first three summands give 



+ 



+ 



Performing the change of variables u = \JE cosh t in the fourth summand and using ( 71 ) , 



Thus, 

yi{Q)dq 



Er= o(i + i)[-JV2]^^' 

24\/n2 - E 
5ET=objA^^' bn A 



du 



12 



^{j + l)bj^2A^ + o{E'). 



+ ■ 



12(^2 _^)3/2 6^(^2_^)l/2 12(^2 _£;)!/ 

+ — + -y6M^'^'+o(i?°). 

12 6S 2 ^ ^ ^ ^ 



i=o 



Note that the error term o{E^) in the previous formula cannot be simply replaced by 0{E^), as terms 
of order ElnE can also be present. 

4.3 Formal monodromy along a[, k even 



The calculation is analogous to the one performed in Sec 4.2 therefore, only the answers will be given. 
For definiteness let us work with k = 2. 



In this section 4.3 denote: 

^ 00 
27ris,>{E,h) = -{AS){E) + Y,h^Ay,){E), 



j=0 



ASiE) = z f ^E - {f'Ydq', Ay,{E) = j y,[q)dq. 



In this section 4.3 we write A = —f'{q2 — e), Oj = 0^(92)) bj = bj{q2)- 
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4.3.1 The AS* summand in 2iTis„i . 

^2 



4.3.2 The Ayo summand in 2TTis^i^ 

We have 



Ai/o = / yo{q)dq = -arccosh 



:b 2 



where the branch of arccosh is chosen so as to coincide with the principal real value of of arccosh for 
-E > and ^2 — £ on the real axis immediately to the left of . 

4.3.3 The Ayi summand in 2Tris„i^. 



We have 



'2 j=0 

4.4 Formal monodromy along 7^ 

The closed paths jk, 7^ were defined on figjij 

Once we calculate s-y^, we will automatically know sy^ since by |Glll Lemma 5.2] 

The following notation will be used in the rest of the paper: 

u:^,{E)= I ^E-{f'{q)fdq. 
J Ik 

Proposition 4.2 We have the following equality of E- dependent formal Laurent series in h: 

00 

27ris,,= '-u^,,{E) + Y,^(^']iE)h^ 

j=0 

where for E ^ 

''^^~~n^) — 4 — 

<^ = -2ni, (25) 
J7« = 0{E), (26) 
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and where 



02 



as in (15). 



Proof. In order to have fewer indices, put k = 1. We need to show (24), (25), and (26). 



Formula (24). Performing the substitution u = —f'{q) (so that q = Qi corresponds to ti = —\JE 



and q = q^ corresponds io u = yE) and using (16), we have 

He 



- {f')^dq = 2 / ' ^ VE-u^iao - am + ogu^ - ...)du = ^ ^ (j + 2)!! ^'^^™ 



71 



2^^7rao + 2^E2^a2 + o{E^) 



J even 
p2 



(27) 



Formula (25). The integration 

f yodq 

J -71 



dq 



ff" 



r 



^2{E-{f'Y) 2i^E-{f'Y 
is performed by substitution u = —f'{q) in the integral. 



-2m 



Formula (26). Make a change of variables u = —f'{q), choose a closed clockw ise con tour 7 in 



the n-plane containing inside it points u = ±vE, and proceed analogously to section 

yi{q)dq = 



4.2.3 



71 



^ 8(^2 - E)5/2 8(m2 - ^)3/2 

Er=2(2-j)(j-i)?'.^^-' 



3-8(n2 -^)i/2 



-du 



27r J] 

2fe=j>2 



^(j-2)(j-l)(j-3)!! 



,_i(2fc-l)!! 



3-8 (j-2)!! 12^ (2fc-4)!! 



0{E^) 



Here use the formula ^j^^du = -27r^^/2('^ for even A;. □ 



4.5 Monodromies from — £ to g^+i — £. 



Define Mj, M'- to be the monodromies of the formal WKB solutions along the paths shown on figjTj 
where Mj are taken on the first sheet of the Riemann surface of p{q) and Mj on the second. Note 
that in [UTT] we denoted M,-, Mj by Aj,A'-. 



Lemma 4.3 We have: 



Mj = exp I ^mJ."^^ + ^ h'^mf^ 
I fc=o 
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contour for M contour for M' 

j j 



Figure 7: Integration contours defining Mj and M'. 



where Wl - = Wl - (E), k > —1, are analytic functions of E at least for < E « 1, such that when 
E-^0, 



{-ly-^iri + - In 



1 



[f{qj+i-e)-f{q,-e)]-E 



f'{qj+i-e) 



E 



E 



4 [f'{q,+i-e)]^-E 8[/'(g,+i - e)]^ 8[f'{q,-e)]^ 



+0{E^ 



Ai) ^ f"{qj+i-e) f"{qj-e) 

2{nq,+i-e)Y 2(/'(g,-e))2+^^ ^' 



where J^^''^^^^^^' means that the integration path lies on the first sheet of the Riemann surface of p{q) 
and In stands for the branch of logarithm that is real for positive arguments. 



Proof. The statement about 



t(-i) 



is obvious since 



VE-[fiqWdq 



92 -£ 



1 E 



<?!-£ 



2 ifiq)) 



:)dq + OiE^). 



We will show the statement about M^*^^ for j = 1- In the integral 



4°) 



(g2-E)i 



(gi-e)/ 



yo{q)dq 



(<?2-e)/ 



rr 



r 



2{E-{f'Y) 2i^E-{f'f 



dq 



the first summand yields 



{l2~e)i 



rr 



In 



to-.), '2{E-U'?) 

E-[f{q2-e)? 
E-[f'{q,-eW 



2TTi 



K E-[f'{q2-er ^ 
4 E-[f'{q,-eW 

ni 1 fE-[fiq2-e)l 
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and the second summand 

cq2-e jii 



qi-e 



92 -e jii 



1 + 



1 E 



E 



dq + 0{E'^ 



2 \f'{qi-e)J 8[/'(g)]2 



q=q2-s 



E 



q=qi-e 



+ 



+ 0{E^ 



2 [-/'(gi-e)] 2 8[/'fe - ' 8[/'((?i - 
If j is even, the passage from complex Ln to the standard branch of logarithm will be slightly different. 



As for 



obtain: 



again for j = 1, we use (21), replace each summand by its limit for E ^ 0, and 



{q=q2-e)i 



{q=qi-£)l 



yi{q)dq 



9=92 -e 
q=qi-£ 



q=q2-e 
q=qi-e 



Hfnr? fu"? 



if 



ll\2 



/(3) 



8(/')5 2U'Y 8(-(/')3) 4(-(/')3) 4(-(/')2) 

"1 y// q=q2~e 



dq+0{E) 



W? 2(/')2 



dq + 0{E) 



2(/ 



i\2 



+ 0{E). 



9=91 -e 



□ 



Analogously, we have: 
Lemma 4.4 14^e have: 



M- = exp 



1 ~ 



(-1) 



h J 



fc=0 



where M 
0, 



(^) - 



(-E), /c > — 1, are analytic functions of E at least for Q < E « 1, such that when 



^ 2 



t(0) 



(9j+i-e)// 



E 



+ 0{E'), 



E 



4 [f'{q^^,-e)f-E ^[f\q,-e)f 8[/'(^,+i - e)]^ 



+ o(ii;2). 



0{E), 
■es on i 

and In stands for the branch of logarithm that is real for positive arguments 



where J^^^^^^^"^^^'' means that the integration path lies on the second sheet of the Riemann surface of 



5 Passage from the equation Pip = Eip to the equation Pip = hErip. 



We saw in Sec 3.3 that as long as we work with formal solutions of = ^"4^ -PV' = hEj.ip, the 
passage from E to hEj- is straightforward. The situation is more subtle once we begin to study the 
correspondence between formal and actual resurgent solution of the equation Pip = Eip or Pip = hErip. 
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This correspondence is expressed in terms of connection formulas and is the main topic of this section 

El 

We will repeat here the formal calculation of the connection coefficients across the double turning 
points of the equation Pip = hErip by the exact matching method, cf. our exposition in [GIU §7] and 
references therein. Compared to |G11] . now we will push the calculation to one more order in h. 

Note that the connection coefficient called cj in [GIU §7] will now be denoted c'j, consistently 
with the notation of [GIU §8]. 



5.1 Stokes curves and Stokes regions for Pip = Eip and Pip = hErip. 



(39) 



The purpose of this subsection 5.1 is to give an extended literature reference for the formulas (30) and 



Recall that following the setup of |G11| we assume < arg/i << 1. In the complex plane of 
q let us draw the Stokes curves - the locus where a discontinuous change of exponential asymptotic 
expansions of solutions of our differential equation may happen. Stokes regions are the domains into 
which the complex plane of q is split by the Stokes curves. 

It is known, e.g. [Y83j , that for a Schrodinger equation 

[-h^d^g+Vo{q) + hViiq)]iP{q,h) = 0, argh = a, \h\ ^ 

with entire Vo{q), Vi{q), the Stokes curves are given by the condition J^^ ^yVo{q')dq' G e*"M, where q^ 
runs over the zeros of Vo{q). 

For the Witten Laplacian, this means the following. In the neighborhood of the real axis in the 
q plane, the Stokes curves for the equation Ptp = Eip, resp., Pip = hErip, look as on fig. [8] and [9) 



real axis 




Figure 8: Stokes curves for the equation Pip = Eip, < E « 1. 




Figure 9: Stokes curves for the equation Pip = Eip, E = 0, or for the equation Pip = hErip, Er £ C. 

As we explained in secj2j we are interested in solving the connection problem for the equation 
Pip = hErip between the intervals {qj-i,qj) and {qj,qjj^i). Looking at the figj9| this is the same as 
solving the connection problem between the Stokes regions containing points qj — e and qj+i — e. 
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For the equation Pip = Eip, < E « 1, the answer to the connection problem between Stokes 
regions containing qj — e and qj+i — e is particularly simple, see |V83] . [DP99] . and our exposition 
in [Gin §7], and will be stated in ([30]) and pOl). 



5.2 Exact matching method around qi. 

We are now going solve the connection problem across a double turning point qj of Pip = hE^ip in 
two representative cases - for the double turning point qi where f{q) has a local minimum in this 
subsection 5.2, and for the double turning point q2 where f{q) has a local maximum in the subsection 



5.3; similar results will hold for other real local extrema of f{q)- 



5.2.1 Definition of the reduced connection coefficient {c'lY {E,h) 



Let be a small positive number. Consider two formal solutions ip+{E,q,h) and ip-{E,q,h) of (11) 
corresponding to the first and to the second sheets of the Riemann surface of p{q) and normalized in 
such a way that ip+{qi — e) = ip-{qi — e) = 1, where e is a small complex number, arge E [0, |], i.e. 



(i?, g, /i) = exp I y - ^E - [r{q')f + y, {q')dq' | (28) 

An actual resurgent function solution (p{E,q,h) of Pip = Eip, < E « 1, whose exponential 
asymptotic representation on the interval (g^ — l,gf ) is given by cp{E,q,h) ~ ip^{E, q, h), has on 
{0.1^12) asymptotic representation 

<P{E,q,h) ~ ilj+{E,q,h) + c'^{E,h)ilj.{E,q,h) (29) 

where 

c[{E,h) = exp(27ris^J, (30) 
cf. figjsjfor the definition of ai and sec 4.2 for the calculation of Sai- 



Since the contour cJi gets pinched by q^ and qf when £^ — t- 0, substituting hEr for E directly 
into the expression of c'i{E) is problematic. Therefore, the following trick is used. 

Introduce the reduced connection coefficient {d^Y'^'^ by 



CliE.h) = ^ {c'^Y''{E,h). (31) 



4.4 



Here s^^ = {E, h) is defined on figj4] and computed in Sec 

It is believed that {c'Y)^'^'^{E,h) can be written in the form exp{^^_-|^ h^A}^{E)^ where A}^{E) 
are holomorphic functions of E near the origin; therefore, {c'lY^'^ {hEj. , h) is straightforward to define. 
Further, s^-^{E,h) is representable by power series in h with coefficients analytic with respect to E, 
therefore also s.y^{hEr, h) = YlV=o ^k{Er)h^ makes sense as a power series in h with coefficients 
analytic functions in E^. Substitution of t = s-y-^{hEr, h) into h^^^ and ^{-t) done by means of the 
Taylor series expansions of the latter two functions at t = BQ{Er). Combining these definitions with 



(31 ) allows us to define c'i{hEr, h). Notice that the expression /i*^2 is representable as an E'j.-dependent 



expansion in powers of h and h In h. 
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Remark 5.1 This passage to the hmit and replacing E > hy hEr has been used in two pa- 
pers |DDP97j and |DP99| . but ideally it would need a more solid mathematical justification. The 
first issue is purely algebraic: one has to show that the coefficients in the asymptotic expansion of 
{c'iY^'^{E,h) with respect to h are analytic functions of E near the origin. This problem is known 
in the literature as the Sato's conjecture. See ^SS06j for its solution for the Schrodinger equation 
with the harmonic oscillator potential. The case of a general potential may perhaps be proven us- 
ing reduction of an arbitrary potential well to a harmonic oscillator using methods of |AKT09j and 
references therein. |^ The other issue is analytical: in our setup, |G1H §2.1], resurgent functions are 
equivalence classes of analytic functions in h and thus the notions of limit, parameter dependence, 
analytic continuation should be transferred from analytic functions to resurgent function with care. 
In particular, the fact that we can define c'i{hEr,h) does not yet formally imply that this quantity 
gives the connection coefficient for the equation Pip = hEri/j- We hope that all these technicalities 
will be resolved in due time. 



5.2.2 Calculation of {c[y'"^{E, h) for < « 1. 

Denote 

oo 

{Ar\E,h) = exp{^ 



(32) 



-1 



In this subsection |5 . 2 . 2| we are going to evaluate 0o, and 0i for — )• 0. 

gives: 

1 



Formula (44) shown is sec 5.4 



27r/i''^i + 2 

r(-s^j 



exp 



27r 27r 



UJ, 



71 

2-K 



-Ln 



UJ. 



71 

2tt 



h- + 0{h^) + 0(Mn E)h } , 



6uj. 



71 



(33) 

where the notation 0{h?) + 0{E\'n. E)h means an ii^-dependent power series YlT=i^k{E)h^ where 
Gk{E) are analytic functions of £; for < £^ << 1 and Gi{E) = 0{Eln E) when E 0+. 



Combining (30), (31), (33), and results of sec 4.2 we have: 



a) 



2-K 2tt 



ijj. 



71 



2-K 



where AS was computed in (20); more explicitly, 

E^ 



e_i = 2[/(gi)-/(gi-e)]-ao^Ln {2A)+aQ^ - a2—E 



1 



j=l or j>3 



A^E+^ — -E'^A^- 
J 4(j - 2) 



^^Ln(2rfe)) + ^ + ^Ln(2r(gi)) " a,^l.n 2A + o{E^) 



b) 



G 







arccosh 



A 



E 



Kl 1^ 

— + -Ln 
2 2 



UJ. 



71 

2k 



■^We thank Shingo Kamimoto for pointing out to us both of these articles. We thank professors Aoki, Kawai, and 
Takei for explaining the result of [AKT09j and its significance. We must confess that as of this writing we have not 
carefully thought through the argument of [SS06| . 
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Ln 2A 



E 
4^2 



vri 1 
~2^2 



c) 



TT 



71 



Ln 



12 



a2f"{qi)E- 



V"{qi)' 



+ o{E). 



1 r(gi-£) , a2[r(</i^^2 

2 (/'(<?! 



+ 



24 



+ o(E°). 



We remark that by lemma 5.6, this last o{E ) is actually 0(E). 



As we pointed out in Remark |5.1[ there must be a conceptual way of proving analyticity of 
@j{E) near £" = for all j > — 1; for now we will prove analyticity of Go, Gi directly in section 
[531 

Note that the infinite sums appearing in the expression for G_i for specific f{q) can be evaluated 
by integration: 



, = 1 -^91 



ao- 



-)dq, 



i=3 



qi-e 



<?i 



f'iQ) 

_J: aof"iQ) _ aifiq) _ a2f"{q) 

[f'{q)f [f'{q)? [f'{<D? f'il) 



dq. 



(34) 
(35) 



5.2.3 Calculation of c[{hEr, h) 

Combining the above formulas, we are now able to to write down an expression for c\{hEr,h) to 
leading orders in h and hlnh valid for any Er £ C. 

Since in the rest of the paper we will be substituting a small resurgent function for Er, i.e. 
Er G C will be replaced by Er{h) — )• for /i — )• 0, we will study the behavior of c'i{hEr,h) for 
Cb Er-^ 0. 

We have the following equality of formal power series in h with coefficients holomorphic functions 

in Er'. 

Er 



s-yi {hEr, h) 



1 



+ 0{Er)h + 0{h' 



where 0{Er)h + 0{h'^) stands for ^l^i Ak{Er)h^ with Ai{Er) = 0{Er) for E,. 0. Therefore, 

E 1 °° 

l^s^,(hEr,h)^ ^ ^-27^-2 e^p[Y^Ak{Er)h'' In h] 

k=l 

and we will write the exponential as 1 + 0{Er)h\nh + 0{h^ In^ h). Here we are obviously facing an 
-Er-dependent expansion in h and hlnh for /i — t- 0. In all our formulas from here on we will treat hlnh 
as dominating h and treat /iln™ h as dominated by 1 for /i — ?■ 0. 

Further, 

_1_ = i + jt + 0{t^), as t^O, 

where 7 = 0.5772... is the Euler-Mascheroni constant, hence the equality of power series in h with 
coefficients analytic functions of E^: 



rf-s. 



^ 00 

= ^BkiEr)h' 



k=0 
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BoiEr) = 1 + 7: 



Er 



-2f"{qi) 

Collecting all our calculations, we have: 
Proposition 5.2 We have 



+ 0{El), Bi{Er) = 0{Er). 



Er 



X exp 



2[/(9i) - fiqi - e)] 



(1 



Er 



f'iqi) 



\ In 



, , h2 1 f"{qi-e) fl2(.r(gi)) 
12 2(/'(gi-£))2 + 24 



2 \ 1 



X ([1 + 0{E'^)] + 0(S^)Wn /i + 0{Er)h + ©(/i^ln ^/i)), 



where aj = aj{qi), 62 = ^2('7i) were defined in section 4-1 



(36) 



Here the error term [1 + 0{Er)\ + 0{Er)h\ii h + 0{Er)h + 0(/i^ln ^/i) stands for an expansion 
Y.k,£>o^kA^r)h''ln '^h with Boo(^r) = 1 + 0(£;2) for Er 0, Bn = 0[Er), Bio = ©(i;^), and for ah 
other nonzero B^^^ either A; > 3 or /c = 2 and £ = 0, 1, 2. 



5.3 Exact matching method around ^2- 



Let < << 1. Consider the following two formal solutions ^+(£', g, /i) and il)-{E,q,h) of (11) 
corresponding to the first and to the second sheets of the Riemann surface of p{q) 

q 00 

^±{E,q,h)=exp{ / y^E-[f'{q')]^ + ^y,iq')dq'} (37) 
Jao-s " ■ n 



where e is a small complex number, arge G [0, 

An actual resurgent function solution (j)(E,q,h) of Ftp = Eip, < E « 1, whose exponential 
asymptotic representation on the interval {qf , g^) is given by (/){E, q, h) ~ iIj-{E, q, h), has on (gj, g^) 
an asymptotic representation 

^{E,q,h) ~ C2{E,h)i^+{E,q,h)+il;_{E,q,h) (38) 

where 

C2{E,h) = exp(27riSo-^), (39) 
cf. figjs] for the definition of o"2 and sec 



4.3 



for the calculation of s„' ■ 

'^2 



Introduce ^^'^{E, h) by the equation 



C2{E,h) = V27r-- -A'''{E,h). (40) 
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It is believed that C2 {E, h) analytically depends on E in the sense that 



ci'\E,h) = ei,Y>{Y.yQ,m 

i=-i 

and Qj{E) are analytic functions of E in the neighborhood of the origin. 



The Stirling formula (44) and (23) give: 



27r/i'^2"^^ 



exp 



1 [ojy^ ^'^^hn ^^"'^ 



2tt 27r 



2tt 



2 V 27r 



TT 



6w. 



+ 0{Eln E)h + 0{h^ 
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(41) 



Similarly to the previous subsection, with A = —f'{q2 — e), aj = aj{q2), bj = 6^(52), we have: 

G_i = z / p{q,E)dq - ^^ + ^^Ln ^"^^^ 
Ja' 



2-K 27r 



27r 



^2 / Eao E^a2 



-2[/fe) - /fe - e)] + i^aoLn (-2^) - + , ^ 



+ 



Ln 



ao 



A^E 



-a2- 



+ ^2 



p2 p2 /I 1 

a2— Ln (-2^) + E ( 7^' + ^, ^''''E ] + 



16 



J 



4(i - 2) 



eo = /,.o(.,^M.-^Ln (IJ) = 
{-2A) , , a2|rfe)|-i? 



-Ln i — Ln 



01 



Ql = [ yi{q,E)dq + ^+o(E' 



+ o(ii;); 



+ 0(^°). 



Using 



we can write 



s^i^{hEr,h) 



Er 



2f"{q2) 



1 



Er 



+ 71 



E?. 



T{-syJ 2f"{q2) ^[2f"{q2)] 



+ hO{Er)+0{h^), 



+ 0{El) + hx 0{Er) + 0{h^). 



We combine the results as follows: 



24 



Proposition 5.3 We have 



X exp 



-/(92) + /(92 - e) 



-(1 



-)ln 



ni2)' V-V"il2) 



=1 J 



(I 1 r(g2-g) «2[r(g2)]^ 

X ([l + 0(£:2)] +0(^^)/iln;x + 0(£;^)/i + 0(/i2in2/j))^ 
where aj = aj{q2), 62 = &2(<?2) w^ere defined in section 4 



(42) 



/(?) 



Supplement. Notice that if Er = 0, then e '» is the formal solution of (11) corresponding 



/(9) 



to the second sheet determination of p(g). But e ^ is also the actual solution of (11), and therefore 
the connection coefficient vanishes: 

C2(0,/l) = 0. 



As a side note, the Supplement allows to collect the error estimate of Prof 5.3 similarly to how 
we did it in Prop j5.2[ 



5.4 An application of the Stirling formula 



In this subsection |5.4| we will derive the formula (44) which we used in (33) and (41). 



Let us first work with a fixed < £' << 1. We have calculated earlier that s. 



71 



2TTh 



where /3i = (3i^i + 0{h) is a power series expansion in h. and /3i,i = 0{E). Since —s-y-^ has a positive 
real part which goes to infinity as /i — t- 0+, we can apply the Stirling formula to r(— s^^) to get 



27r/i''^i + 2 



A few routine steps of simplification bring us to 



27r-(l-3^S7/ + 0(/i2)) 



/2^/l'^l + 2 



exp{-^} 



r(-s. 



71/ 



2-K ) 



nh 



-{l + 0{h')). 



) 



(43) 



In (43), the LHS is a true function, and the RHS its asymptotic expansion valid for E > and for h 
in a small sectorial neighborhood of zero in the positive real direction. 

Now let us remember the -E-dependence in (43). Then (1 + 0{h^)) on the RHS should be 
interpreted as a power series 1 + '}2'k^=2^k{E)h^ where Gk{E) are analytic functions of E for < 
E « 1. 



From sec 



4.4 



wejcnow that /3i = j3i{E,h) = J2'k=o^k{E)h'' with Mo{E) = 0{E). We can drop 

m 

l + 0{h^) + 0{h)0{E\iiE)), (44) 



/?! from the RHS of (43) at the expense of changing the error term as follows: 

exp{-^} 



27r/i^Ti + 2 



'71/ 



27r . 



"71 I 1 
■ 2TTh """2 Q 
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where 1 + O^h?) + 0{h)0{Eln E) stands for an £'-dependent power series 1 + YlV=i '^'ki^)^^ where 
are analytic functions of for < ^ << 1 and G'^{E) = 0{E\n E) when E ^ 0+. 

5.5 Partial proof of analyticity of {c'i)red 



The purpose of this subsection 5.5 is not to give a full justification of the methods, but rather to prove 



some results confirming that our approach is consistent and makes sense. We are going to show that 



G_i(i?), Q)q{E), Qi{E) from (32) are indeed analytic functions of E. 
Lemma 5.4 The quantity 

= A5+^-^ln[-^] 
27r 27r ^ 27r^ 

is analytic with respect to E in the neighborhood of zero. 

Proof. We need to show that the term containing In E in A5 cancels -^Ln E. Indeed, we 
have seen that 

oo ..^ oo „o 

= -2^{-lfak j vf'^/v? - Edu = -2'^{-lfak / E 2+^ cosh^ t sinh^ tdt 

J A j^_Q J arccosh A/\/E 

00 „o 

= _2^(-l)'=afc / 2-^-2i?t+i(e* + e-*)^(e* - e-'fdt. 

./arccosh A/y/E 

Writing the integrand as the sum of exponents and integrating, we realize that only the summand 
e^^dt will eventually give rise to a logarithmic singularity for £' — >• 0. Writing reg{E) for an arbitrary 
function that is analytic with respect to E near the origin, we have: 

00 „o 

e_i = regiE) - 2 ^ a^, / 2-'^-'E^+Hcij' - 2Ci^ + Cil')dt 

k=2j=0 -/arccosh A/\/E 



reg{E) - a2j2-^^ E^+^^{1 - ^—)arccosh { 



k=2j=0 

reg 



k=2j=0 ^ ' 



The singularity that comes out of E is the same by formula (27). □ 

Lemma 5.5 Gq is analytic for E around 0. 

□ 



Proof is obvious from (70) and proposition 4.2 



Lemma 5.6 0i is analytic for E around 0. 
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Proof. The fact that Gi has no pole (i.e. singularity has been demonstrated in section 5.2 
Now let us check that all logarithmic singularities E^'hn E are absent in 0i as well. The question 
reduced to identifying the logarithmic singularity in the integral along the contour a, figj6j 

^ Ef=o(j + l)[-jV2]n^' ^ 

= Tl^E^^(^' + l)H^i+2] / COSy(i)dt = 

~^ J arccosh iA/%/E) 

regiE) + ^ ^ E^j + l)[-jbj+2]2-^ C^'^^\-aiccosh (A/Ve)) 



j=2k>0 



The logarithmic singularity in the term of Ln w^^^ ^ ^ comes from — /3Ln { — '-^), i-e. from 



/3Ln E, where 

fc^i (2fc-l)!! 

27ry,/^^^ 24^"^^" (2A;-4)!! 



1 /■ -j^ ^ 

That means that ii^'^Ln E terms in ©i cancel for all k. □ 



5.6 Remarks on more precise calculation of c'i{hEr,h), C2{hEr,h). 



Let us discuss what needs to be done to calculate more terms of c'i{hEr,h), both with respect to h 
and with respect to E^; similar ideas apply to C2{hEr, h). 

In order to determine c'{hEr,h), we had to perform three pieces of somewhat laborious cal- 
culations - computing Saj^{E,h), S'y^{E^h), and the expansion 43 all other steps were more or less 



substitution and cancelation of terms. An analog of (43) is available by the Stirling formula to any 
order in h. Computations of Sai{E,h) and s^^{E,h) are parallel to each other, so let us discuss the 
former. 



of the type 



/(H)(g).../(^f)(g) 



The j-th term yj of the series (|12[) solving the equation ( 13 ) is a linear combination of fractions 



where f^'"'\q) denotes the ip-ih derivative of f{q), and A; E N. 



[£;-(/'(g))2]'=/2 

Similarly to (17), we can write f^^^^q) as a power series in f'{q). Thus, the integral Ayj, cf. 

f'{q) into an integral of the type ^^^^ " 



(19), can be transformed by substitution u 



(£;_„2)fc/2 



where a is the contour on figjGj If k is even, this integral can be taken by the residue formula. If 
k is odd, this integral can be calculated using several steps of integration by part s to r educe the 
denominator to \/u'^ — E and then using substitution u = ^/E cosht, similarly to sec 



4.2.3 



In terms 



of the coefficients Cm, the formula for Ayj{E) can be written precisely as a function of E. 
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6 Ingredients of the quantization condition for the equation P?/' = 
hEri/j. 



The purpose of this section [6] is to calculate the ingredients of the quantization condition mentioned 



in Step 3 of sec 2.1 



6.1 Calculation of fj.j{hEr,h) 

Define: 

/i-i -1 2iTis^ ■ (E,h) -11 II —2iTis-y .(E,h) 

(j,j{h,h) = 1 — e , J odd; fij = 1 + e ^ ,J even 

for E G C, \E\ sufficiently small; ^j{hEr, h) is then defined by substitution of hE^ in place of E. 
Using Prop |4^ and performing routine simplifications, we obtain 

^J.,{hEr,h) = -^(1 - '''^^^^('^''^\ Er - ^j^frW + 0{E?.)] + 0{h')), J odd, (45) 



\n,,)r' 4 hEr-^^^E.m + OiE^^)] + Oih^)) je.en (46) 

(in the sense that the term of the power series with respect to h is given up to 0{E'^) for E^ — )• 0, 
term is given precisely, and h? and smaller terms are neglected). 

6.2 Calculation of rj (/ii?r) ^) 
We define for ^ G C 

for odd j: Tj = Tj{E,h) = CjCj+iM~^Mj, 
for even J : tj = Tj{E,h) = CjCj_^_iMj{Mj)~^ , 



and the definition is then extendable to Tj(hEr, h) by means of Props, 5.2 and 5.3 
Let us now evaluate two representative cases of Tj 



Calculation of ri. Using lemmas 4.3, 4.4 and propositions 5.2, 5.3, and inserting E = hEr 
into the corresponding formulas, obtain after routine simplification: 



,f(n)-f(i2) 



ErIT 



Er 



Er 



X exp 



.7 = 1 ^ ^ 



-In 



/"(ft) vm^) 



+ 



X exp 



X exp 



Er 



{Q2-e)i 



+ 



TTlEr 



12 24 
^2(92) a2{q2)[nq2)? 



12 



24 



f'{q) f"{qi) 



([1 + 0{E^r)] + 0{Er)h\n h + 0{Er)h + 0{h' In^ h)). 
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Here fl^g^_^-ll denotes the integral along a path lying within the domain of definition of figjsj 
Formulas for other tj with odd j are analogous. 



Er 



Er 



Er 



X exp 



exp 



Calculation of T2. Analogously, 

T2(hEr,h) = e2^^^^¥^ ^ h^rW^Ul+^ - +^ - ) 

a2iq2)[f"iq2)? 



i=i 



24 



X exp < —Er 



TViEr 



12 24 



and similarly for other tj with even j. 



([1 + 0(^;)] + 0{Er)h\n h + 0(£;^)/i + 0{h^ In^ /i)), 



In a calculation of these monodromies for a specific f{q) we can use formula (34). 



6.3 Formula for the product of Mjs 

Let us calculate the expression that we called 1 + ErK{Er, h) in |G1H Lemma 8.2]: 
l + ErK{Er,h) =^ M[{hEr,h)M^{hEr,h)...M^^{hEr,h) = 



(by lemma 4.4 ) 

V ^ J2nd sheet / WJ/ 

where ^gnd sheet ™6ans that the integral is taken along a loop from some go to go + 1 along a loop lying 
within the domain of definition of 6-. 



7 Quantization condition and formulas for eigenfunctions. 

In this section we give explicit details and formulas to execute the method of secj2] of calculating 
low-lying eigenvalues and corresponding eigenfunctions of the Witten Laplacian 



7.1 Quantization condition 

The explicit form of the transfer matrix F{Er, h) from Q was discussed in jGlH §8]; the quantization 
condition for the Witten Laplacian will be written in terms of a related 2x2 matrix Go = Go (-En h) 
whose coefficients are E^^-dependent resurgent symbols in /i, or, by abuse of notation, Ej.-dependent 
exponential asymptotic expansions in h of the form 

j 
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with Aj^k,£{Er) analytic functions of Er and cj G M independent of Er. 
The quantization condition is: 

1 



ker(Go 



1 + ErK{Er, h) 



Id) ^{0}, 



or equivalently 



1 

1 + ErU 



+ Tr Go - (1 + ErK) det Go = 0, 



(45 



(49) 



where the coefficient K{Er,h) was defined and computed in (47). In the case when f{q) has n local 



minima and n local maxima on the period q G [0,1), the matrix Go = Go{Er,h) has the following 
explicit form: 



Go 



n 



k=n,...A 



where the factors in the product are arranged, from left to right, in the order k = n,n — 1, 1. In 
particular, if n = 2, 

_ f T4{t^^ + l)T2{r^^ + 1) + ^4(^3^3"^ + l)(M2Tr^ + 1) 'i~4{t^^ + l)r2(^irj"^ + 1) + + 1)(miM2Ti~^ + 1) 

(50) 



where fij = ^j{hEr,h) and tj = Tj{hEr,h) were calculated in sec 6.1 and 6.2, respectively. 



7.2 Formulas for the eigenfunctions 



Assume Er{h) is an exponentially small solution of the quantization condition (49) 



Let (j)±{hEr{h), q,h) be formal WKB solutions obtained from (j)±{hEr,q,h) of sec 3.2 by sub- 
stituting Er{h) in place oi Er & C into the expression of (p±{hEr,q,h). As Er{h) has an exponential 

asymptotic expansion with several, usually countably many term, Er{h) ~ aj{h,hlnh), the 

g-dependent exponential asymptotic expansions of (j)±{hEr{h),q, h) will be of the same nature. 

Denote BoiE,h) = 



Let go be as in sec 



3.2 



= ME, qi - e, h), B'oiE, h) = <P^{E, qi - e, h). 

Suppose a vector of resurg ent symbols (z|°\ zi°^)^ belongs to the kernel Then, cf. fGTTl 

^8.2], the vector 



D 
D 



(0) 

+ 

(0) 



i?^i?o-i(c;)-izf 



z 



(0) 



belongs to ker {F{Er, h) — Id) = and thus the eigenfunction corresponding to the resurgent eigenvalue 
hEr{h) will be represent able, for q G (0,gi), by an exponential asymptotic expansion 



{h)(t)+{hEr{h),q,h) + D"-^' {h)^-{hEr{h),q,h) 



(0), 



The connection formulas for the equation = hErip can be rephrased as follows. For q G 
{qj,qj+i), j = l,...,2n, g2n+i := Qi + 1, the same eigenfunction is representable by the following 
exponential asymptotic expansion 



of {h)<^+{hEr{h),q,h) + D^l\h)(j)^{hEr{h),q,h), 
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Df{h) 



and are defined recursively based on Z. 



B'^ihErih), h) ^-yt M[{hEr{h), h) r (c'AhErih), h))-\ j odd 



Bo{hEr{h), h) n Me{hEr{h), h) 



Cj{hEr{h), h) 



J even 



Z^\h), 



(0). 



z^'^ + ^.z^'-'^ 



j odd; 



z 



z 



+ 



r7\z^r'^ + Z^r'^ 



J even. 
(51) 

The same formulas will look nicer if we introduce the following basis of formal WKB solutions. 
Let tp-{E, q, h) = cj)^{E, q, h), in particular, tp-{E, h, go) = 1- Suppose is in a neighborhood of the 
positive real axis; accordingly to our conventions, we think of 'tp-{E,h,q) as defined on the second 
sheet of the Riemann surface of p{q)- Then ip-{E,q,h) can be analytically continued (as a formal 
WKB solution) across the cut connecting qi{E) and qfiE) along the path a[, figjsj we thus obtain a 
formal WKB solution tp+{E, q, h) on the first sheet of p[q). After that, substitute hE,., where E.^ G C, 
or hEr{h), where Er{h) is a small resurgent function, to obtain ip±(hEr,q,h), resp. ^±{hEr{h),q,h). 

We can write the exponential asymptotics of the eigenfunction corresponding to the eigenvalue 
hEr{h) as 



Df{h),lj+{hEr{h),q,h) + D'l\h)ilj.{hEr{h),q,h) on ((?j,g,+i), j = l,...,2n. 



where 



D 



Z^ X 



n 

£=1 



.(-1/ 



D 



z 



(52) 



, /(g)-2/(9i)+/(gro) 

~e h 



Remark that 

ijj^{hEr, q, h)'^e ^ h ~^ Tjj^{hEr,q,h) 

in the sense of the notation ~ introduced on p|3j 

We also see that removing the condition ip-{E,qo,h) = 1 will only rescale both solutions 
'^±{E, qo, h). Thus, up to rescaling the vector [il)j^{E,q,h)^il)^{E^q,h))^ we can make our formulas 
independent of the choice of go- 



Remark also that if E^ satisfies the quantization condition, then the corresponding eigenfunction 
of the Witten Laplacian is periodic and we must have 



[MiM2M3M4]-^D. 



in(0) 



of'; [M[M'2M^^M'^]-^D 



' Ti/r' ]\/r' A/f'l-l nC^) 



D 



(4) 



(7) 

Rewriting this condition in terms of Z^ , we arrive at 



Z 



(0) 



l + ErK)c[c4{Miy^M^z':^^; Z^^> = {l + ErK)Z'L>. 



(0) 



.(4) 



(53) 



(7) 

If Z^ are calculated without algebraic mistakes, they must satisfy the formulas (53). 
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8 Example 1. 



Let 



fiQ) 



1 

2^ 



sin27r(q + -) + cos47r((7 + -) 



(54) 



/'(g) = cos2^(g+^)-2sin4^(g+^; 

The summand g in the argument assures that /'(O) 7^ 0. 
The critical points of f{q) are: 



qi = I 

92 = i - ^ arcsin ^ 

„ 5 

93-8 

94 = i + ^ arcsin \ 



fiQi) 

f{Q2) 



= 
9 

TV 

9 
IGtt 



f'iqi) = 6vr, 
r((?2) = -7.57r, 

r(«3) = lOvr, 
r((?4) = -7.5^. 



Now we are going to exploit the symmetry of our choice of f{q). 

Lemma 8.1 Suppose f{q) has two local minima qi, q^ and two local maxima (72)94 o^i^-d satisfies 
f{q)=f{2qz-q). Then 

Proof, a) Observe that for E real, the equation Pip = Etp has two real solutions, i.e. those 
satisfying 

^(9) = W), (55) 

therefore the same equality must be satisfied by any solution of this equation. Thus: 
reflecting a contour 5 on the Riemann surface of p{q) with respect to the (preimage of the) real axis 
of the complex plane of q, while keeping the contour on the same sheet of the Riemann surface of p(q) 
changes the monodromy of a formal solution by complex conjugation. 



Indeed, the monodromy changes from ^l^lgjj 
oiS. 



to 



where 6{0) and (5(1) are the endpoints 




Figure 10: Contours formal monodromies along which equal to Tj{E, h) for < << 1. 

h) The coefficients 'Tj{E^ h) can be equivalently defined as monodromies of formal solutions along 



contours shown on fig 10 Let us show that for > and h real, Tj{E,h) is real. E.g., for T2{E,h) 
perform the transformation of the contour shown on fig 11 : the first transformation is done using |G1H 
Lemma 5.4], the second simply reverses the direction of the contour, and the third reflects the contour 
with respect to the real axis, cf. part a) of this proof. Since we returned to the same contour, we 
conclude that T2{E,h) = T2{E,h). 
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Figure 11: Formal monodromy T2 = T2{E,h) and the related monodromies - step b) in the proof of 
lemma 18.11 

c) The equation P^jJ = Etp is preserved if we replace q by 2q^ — q, and so are the determinations 
of p{q) on the sheets of the Riemann surface of p{q). Hence, a formal monodromy along a contour 6 
coincides with a formal monodromy along a contour obtained from (5 by g 1— t- 2^3 — q. But the contour 
for T3 is obtained from the contour for T2 by means oi q 2^3 — q followed by reflection with respect 
to the real axis. Using steps a) and b), we conclude that T3 = T2- 

d) Similarly, s^^ = s-y^ and hence fJ-2 = fJ-i- ^ 



Solving (49) for Er begins by writing out the TrGo term. Using (50), we have: 
TrGo = 



.-1 



T4T^^T2T{'^ 
+^1^4T-2T3"Vf ^ +lliT2T^^ +/UlT2rf 



+T4T2T3 +TaT2T^ 



-1 



+/ilAi2/t^3/^4T; 



.-1 



+ T2Ti 

+T2 
+ 1. 



(56) 



Using the formulas for r^, we find 

Ti = T4~e~ 

and therefore we can write, loosely, 



9 ^ 



TrG, 







+Ere8^h 

_9_ 

+Ere&^h 
+E^e8'^h 



T2 = T3~e a^hEr, 



1 +Ere 



9 

" S-rrh 



+Ere 



25 



^ 9 



■n 16 

■n 25 

+Eres^h 



+E^e 
+Ere 
+Ere 



34 

9 

25 
+ 1, 



(57) 



by which we mean that, e.g., the exponential type of the summand t/iT2T^ ^ is the same as the expo- 

9 

nential type of Ere~s^ and that, therefore, this summan d co ntributes to the points corresponding to 



E^e svrh^ /c > 1, in the Newton polygon, cf. step 4 of sec 2.1, of the quantization condition (49). 



By |(j1H §8], the first and the third summands on the LHS of (49) contribute only terms of 
type E^, k >1 (degree zero with respect to eh). 
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degree 

ine^^ 

34 
8tt 



25 
8tt 



16 
8tt 

9 
8tt 



1,25 



1,16 



1,9 



a 


a 


a V 


2,34 


3,34 


4,34 


a 


a 


a 


2,25 


3,25 


4,25 


a 


a 


a 


- 2,16 


, 3,16 


4,16 




a 




.^2,9 


. 3,9 


,^9 



degree in E 



Figure 12: Newton polygon of the quantization condition (49) for the Example 1. 
8.1 Solving the quantization condition 

Since f{q) has two real local minima and two local maxima on its period, the Witten Laplacian will 
have two resurgent exponentially small eigenvalues, one of them equal to zero; from now on hEr{h) will 
denote the other, nonzero exponentially small eigenvalue. In this subsection we are going to calculate 
the beginning of the exponential asymptotic expansion of Er{h) using the Newton polygon, cf. step 4 



of secj2l]and [GTTl §9]. 



The upper left portion of the Newton polygon of the quantization condition ( 49 ) for our example 



is shown on fig 12 It is convenient to modify the notation of (10) by putting ^ under the exponent; 
let us write 

J^{Er,h):=^E'^es^ak,c{h) = 0. (58) 

fc,c 



In |G11| we explained that the leading exponential summand of the exponential asymptotic 
expansion of Er{h) is obtained by looking the the north-west edges of the Newton polygon, in our 
case that means - by solving for Eq = Eo{h) the equation 

25 2 34 

then Er will be equal to Eq up to terms of smaller exponential type with respect to h, i.e., Er{h) ~ 
~ a2'34W ^~^' next term in the exponential asymptotic expansion of Er{h), let us make 

a substitution 



Er{h) 



-^ + E,{h)]e'^^ 

^2,34 



and solve the quantization condition for Ei{h) under additional requirement that Ei{h) should be 
exponentially small. 



34 



b 
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16 
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^,-16 






\ 




Figure 13: Newton polygon for Ei. 

In other words, if the quantization condition as an equation on Er{h) was written in the form 
(58), we rewrite it in the form 



Ti{Ei,h) :=Y,E'les^bk,c{h) = 



(59) 



k,c 



where 



J^iiEi, h) = e'^Fi-"^ + El, h). 

02,34 



Let us pfot the (possibly) nonzero summands of the RHS of ( 59 ) on fig, 13 there are no nonzero 
terms to the right of the slanted dotted line through 62 and 63-9. 



We have: 



^1,0 — — «i,25) ^2,0 — 02,34, 

, 01,2501,16 , .0125^2 /'^1,25n3 

00,-9 = h 02,25( ) - a3,34( ) • 

02,34 02,34 02,34 



Analogously to the first step, from the Newton polygon on fig 13 we infer that 



El = e~s^ — ^ + E2 
V bifl 

where E2 has to be exponentially small. The next step of this procedure and a Newton polygon for 
E2 (which we will not draw here) yields E2 G . 

Let calculate the coefficients 02,34, 03,34, 01,25, 02,25, and oi^ie- Only the following four summands 
in Tr Gq contribute to these coefficients: 



34 



E^TT^ 



V/"(9i)-irfe)i-/"(93)-i/"(g4)i 



xh 2 yf'iii)'^ \f"(i2)r f'M^ \f"(u)\J {[1 + 0{Er)] + Oihlnh)); 
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25 



ErTT 



Vf"iQ3)\f"iQ2)\ 



Er 1 I L_ 



/i 2 Vl/"te)r/"{93)y([l + 0(Sr)] +0(/iln/i)); 



25 



-1 ^ 



h''^^i"^(''3)^\i"^^'iA)\) {[l + 0{Er.)]+0{h\nh))- 



Er f 1 1 

— ill ir.„\~TrT 



V/"(9l)/"('/3) 

where the notation (1 + 0{Er)) + 0(/i In /i) means an expression of the form ^Mki{Er)h''ln'^ h with 
Boo(-E'r) = 1 + 0{Er) for Er ^ ^ and ah other = if /i^ In^ /i dominates /iln/i. 

We have: 

^2 1 



TT 



12,34 



V/"(9i)-l/"fe)|-/"fe)-|/"(g4)| 



15V15 



^3,34 — 



V/"(5i)-l/"(92)l-/"(g3)-|/"(g4)| 



15V15 



1 



1 _^ 2 ^ 
Gtt 7.57r IOtt 
vr 



1111 

f"{qi) ^ l/"(g2)| ^ /"(gs) ^ |/"(g4)| 
In ft + 0(/i") _ 4 



In h + 0{h'^) 
2 



152V157r 



In h + 0{h°); 



IT 



01,25 



02,25 



1 1 



+ Oih) 



+ 0{h); 

In h + 0{h'^) 



v//"(a3)|/"fe)| \ f"ii3) l/"(g2)iy Vf"iiB)\f"M\ V/"(93) l/''^) 

21 2 In h + 0(/i") 
"^^10^ 15^^ 2 



■ —In ft + 0(/i°) = ^lnft + 0(ft°); 



^1,16 



+ 0{h) 



0{h). 



We think that the imaginary value of 01,15 is attributable to our choice of arg /i > 0. Passage to the 
case arg/i = would involve techniques studied, e.g., in [M99] . 

Proceed with the calculation: 



r :- 



«1,25 



02,34 



-(.VfiQi) ■ l/"(94)| + Vf'iQi) ■ l/"fe)|) + 0{h) = 6V5 + 0{h). 



60,-9 = 01,16?" + 02,25?" + 03,34r 
4 



18^/3 
Svr 



In /i + 0(/i°); 



-•1,0 



"^1,25 



0{h). 



Collecting the results, we find the following formula for the nonzero low-lying eigenvalue: 



hEr{h) = he 



/ 01 25 



he siTh 



9- / 01,25 



03,24 



+ e 



02,34 

bo-9 
bi,o 



+ Ei 



e-8^(6\/5/i + o(^^)) + e-8^^( hln h + 0{h^)) + £' 

TT 



25 



(60) 
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8.2 Asymptotic expansion of the eigenfunction corresponding to the nonzero low- 
lying eigenvalue. 



Let us calculate the eigenfunction corresponding to the eigenvalue (60). 

As explained in secjrj we will start with {Z^\ Z^)'^ S Ker{Go — j^rj^Id). We can take 

4°^ = A-[Go]i2, zi°) = -A.([Go]ii-— ^) 

where A = A(/i) is any nonzero resurgent symbol which we will choose as A = ^f^^ — to simplify the 

T3 T4/13 

formulas. 



Usmg (57), 

^ 9_ ^ __25_ ^ 34 

Ti = T4~e sirh Ey'^e ; r2 = T3 ~ e sirh Er'^e stJi . 

Let us write down the exponential orders of the various summands in [Go]ii — (1 + Er-n)'^ and 
in [Go] 12- Namely, 

9_ _ . 34 

/^e Sirh rv^e Stt/i rv^e Sivh 52 



[Go]ll - (1 + ^rK) ^ = (r4r2Ti V3 ^ - (1 + Sr/c) ^) + ^4X3X3 ^ + r4T2ri 1 + %^ 



- TjT^ V3 V3/i2 + r4r^ V2 + T3 ^4/^3 + , 

16 9 7 . - 18 



hence 



25 



[Go]ii-(l + i?,K)~i _i ^ r4Tf V2 + (r4r2rf ^3-^ - (1 + Erk)-'] 
n = /^2+i+ ^1 ] +0 8^. 

7"3 r4/i3 ^-v— ' V 7"3 r4^3 



9 



16 

:.e Sirh 



Further, 



9_ ^ 



[Go]i2 = MiT^ r2r3 T4 + T2r3 r4 + /iiT^ r2r4+ r2r4 

+ Ml/"2/»3Tr^^3~^7-4 + /^3T3" V4 + /MJJ^^^ + T4 . 

7 7 18 . - 18 



hence 



[Go]l2 _ 1 , -1 , /^IT"! '^2 



25 



- = 1 + /il^2Tl H ^ \-£ 8^h. 

/^3T3'r4 ' ' ^ ^3 



16 



Thus, 
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(i) 

Before writing down the explicit expressions for , let us derive the following consequence of 



the quantization condition. Using the explicit form of TV Go, assuming satisfies (49), and keeping 



only the largest terms in ( 49 ) , we obtain 



which, taking into account ri = T4 and /X2 = /i4, simplifies to 



-1 n--^ 



(61) 



Now, applying formulas (51), we obtain successively: 



.(1) 
.(1) 



T^l ^2 + ^1^2^^! H ^-1 ^ ^ 



1 - /il + h f 8- 



-^1 V2 



-1 + 



(1 +f-8^ 



(2) 

+ 

(2) 



Z 



(1) 



n +^1 ^2 



-1 + 



16 



16 
Stt 



(2) 

Remark 8.2 It is interesting to note that in the calculation of Z_ the contributions from the leading 
exponential orders in Z^^ and Z^'^ cancel and the nonzero value of Z^^ is due purely to subdomi- 
nant exponentials in Z^^\ Neglecting these subdominant terms would make the rest of the calculation 
impossible. This little algebraic detail is philosophically important: it shows that constructing asymp- 
totic expansions of an eigenfunction of the Witten Laplacian on all intervals {qj,qj+i) must be difficult 
without methods of resurgent analysis. 

Let us continue: 



.(3) 

A3) 



'''' [ rrV2-(l + ^-) 



_9_^ 

Stt 



r2rr^Z« 



-l + £-^ 

-rV2;^+£:° 



where we have used (61) in the last step. Finally. 

.(4) \ 



z\ 
z 



(4) 



)+rrV2;ii+^ 







-T2T^ Z 



/i2T3 +£s^ 
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We finish this section by computing the coefficients D^i from (52): 



To the leading term, 



9 

V ' 

'?_ 



9 



and we obtain the result claimed in ([3]). 



(62) 



9 Example 2. 



Let f{q) be a trigonometric polynomial with two local minima and two local maxima qi,q3 and two 
local maxima 94 on the period [0, 1], where < qi < q2 < q-s < < 1- Up to shifting g by a 
constant we can assume that qi is the global minimum of /, and up to changing f{q) into f{2qi — q), 
that q2 is its global maximum. Changing further f{q) by an affine linear transformation / i— )■ Af + B, 



14 



All these 



we can assume f{qi) = 0, f{q2) = ^, fiqs) = |, fiQi) = f , where < 6 < a < 1, figure 
transformations of / produce easily controllable changes in the eigenvalues and eigenfunctions of the 
Witten Laplacian. 




Figure 14: Graph of f{q) in Example 2. 



We will actually assume that the inequalities are strict: 

Assume: < 6 < a < 1, (63) 
and will gradually put more restrictions on a and b more specific as we progress through this section. 
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Efe " 


1 


E] 


1 














Figure 15: The Newton polygon of (49), situation of Example 2 



In our situation 



1^ l — b a — b a 



In order to find the two low- lying eigenvalues of the Witten Laplacian (one of which equals, as 
we know already, to zero), we need to solve the same quantization condition (49) and will use the 
formula (56) for Tr Gq. 

In the loose sense explained in the Example 1, we have now 



TrGo 



1 +Ere'h +Ere~ +^^e~ 



l-b _ b_ \-a 

+Ere h +Ere h +Ere h 



+Ereh +Ere 



+Ereh 



l + a-b 



+Ef:e—fi 



JL a-b 

+Ereh +Ere h 



+Ere h 

b-l 

+Ere h 
+1. 



l + a-b 



The Newton polygon corresponding to (49 ) is as shown on fig, 15 with the E^e ^ -term coming 
from the /^i^2/^3A^4''"r^'^3~^ summand, and the E^eh term - from the /ii/X2''"r^ summand. We conclude 

h — a 

that the nonzero low-lying eigenvalue will have the exponential type E^'^e^i^ . 

As Go — (1 + -E'r't)"^-^'^ is a 2 X 2 matrix of rank 1, a nonzero vector in its kernel is proportional 
to ([Go - (1 + Erk)-\2, -[Go - (1 + Erk)-\i)'^, i.e. to ([Go]l2, -[Go]ii + (1 + Erk)-^)^. 

For Er'^e~ , the exponential types of various summands in [Go — (1 + Erk)~^]ii, [Go]i2 are as 



follows: 



[Go - (1 + Erk) ^]ii = {t4T2T^ V3 ^ - 1) + T4T2T^ ^ + T4T2T^ ^ + T4T2 + 



£b-a 



-1 + b-a 



-l+3i)-3a 



■ V3 V3Ai2 + T4T1 V2 + T3 ^4/^3 + T4 



l-g 



l+6-2a 



b-2a 



(The first summand should typically be ~e h ^ but it is conceivable that its exponential type is 
actually smaller for a special choice of /.) 



[Go] 



12 
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Figure 16: Values of (a, 6) satisfying, from left to right, (64) 



(65) 



(66) 



(67) 



l+b-2a 



l + 2b~3a 



b — 2a 



The two largest summands in the above formulas are thus ^2/^37"]^ + /i2T^ ^T4 and 

IJ'ifJ-2lJ-3T^^T^^T4 + fiiH2Ti^T4, respectively; so it is reasonable to take 



Z 



(0) 

.to) 



IJ.2T1 V3 V4(/i3 + r3) 



0J12 



[Go]ii + (l + ^rA;)- 



There are too many summands in the entries of Go for us to be able to get an enlightening 
exposition, so we will artificially impose additional assumptions on (a, 6). These assumptions will 
help us select the dominant exponential, the first subdominant, the second subdominant, etc, terms 
in every exponential asymptotic expansion we are going to write down in a moment. There might be 
a combinatorial structure to various inequalities between (a, b) we are going to introduce, but we are 
not ready to comment on it at the present time. 



Under 



additional assumptions: a < 2b; 6 > 2a — 1 



(64) 



(first part of the figure 16), we can write 



[Go]ll-(l + ^r'^)-' 



H2T1 V3 V4(/U3 + T3) 

[G'o]i2 

fJ'2Ti^Tr^^U{fl3 + T3) 

Restricting further to 



1 + 



/i2'ri V4//3T-3 ^ 



+ £ 



-l+2b-a 



1 + 



n 



+■ 



T2 



T"2T3 



/^2Pt3 /^2/X5 



-l+26-a 



additional assumptions: 2a < 36; 5 > 2a — 1, 



(65) 



see the second part of the figure 65, we can absorb the boxed term into the error £ 
We conclude that 



(0) 

+ 

(0) 



'l'2 _ T2T-j, _|_ g 



'1+26-a 



m (^1 + 

]_ _^ {riT2T-'T^'-(l+Erk)'^) ^ g- 
t^2T^^T4ll3T^^ 



l+2fe-aj 
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Using (fsTl), 



Z 



.(1) 
(1) 



T2 _ (r4T2r^ V3 ^-(1+-Erfc) ^) _ r2T3 _|_ 

-1 + ^1 



2f)-a 



with the following exponential orders of the ingredients of Z^^^ : 



T2 . -i+g {t/^T^T^ To - (1 + -Er^) ^) . -l+fc T2r3 ^ -1+6 



e ; 



We have 



(1) 



^-17(1) 



7"2/"l 



hence by (51 ) 



(2) 
+ 
.(2 



^1 / T2 (r4T2ri V3 ^ - (1 + E'rfc) ^) T2TZ 



/^2Aii 



£■-1+26-0 



1 + 

V •^2^1 



gb-2a 



where the coefficient in front of this vector nothing but ^Z^\ 
The formula (51) gives 



7(3) 
7(3) 



Ml ( T2 {T4T2T-i^^T,i^ - (1+ Erk) ^) r2r3 ^^_l + 26-a 



Tl VM2M3 



At2Ti V4/I3T3 ^ 



M2M3 



r2+£:-i+" + M3 [m2-^^] 
r2+£ 



l + 2b— 2a I Ti I c^>— 2a 



Under 



additional assumption: b + a < 1, 2a < 3b 



(66) 



we have T2 G and therefore 



A3) 

7(3) 



A*! T-2 _ (-r4r2ri V3 ^ - (1 + Erk) ^) _ r2r3 ^ ^-1+26-a 



2b-2a 



b-2a 



nM2M3 

7-2/^1 

Tlf 2 



+ r2" +£:^ 



U2 - ^JJ^aia +£:*-2° 

^ T2;ii 



Finally, use (51) to obtain: 



A*! / T2 _ (r4r2r^ V3 ^ - (1 + Erk) ^) _ T2T3 ^ ^-l + 2i,-a 
Tl VM2M3 



M2T1 V4/i3T"3 ^ 



M2M3 



r3-'(M2/i3+f''""")+M2+f' 



M4r3-^(r2 + M3 [m2 - + f^''-^^") + /i2 



T1A'2M3 _|_ £;6-2a 



Under one more 



additional assumption: a < -; 2a < 36, 



(67) 
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we have fJ^^T^ "'^r2~e ^ € <f "^"j and the expression simphfies: 



7(4) 



4) 



|4 



Ml 



T2 _ (r4r2r~V3~' - (l + g,fc)-') 



T2TZ I ^_i + 2b-a 



T1M2M3 



We will see now that the bracket [^t^^ + 1] in the expression for zi^^ is not ~eh as would 
appear from the first glance, but is of a smaller exponential type. Indeed, the quantization condition 
(49) and the explicit form (56) of Tr Go imply 

(/"STa"^ + l)^2T4Tf ^ + /il/U2Tf ^(1 + ^i■■i^^AT^^) = 
(1 + /i3/U4r3-l) = -(/isTs-l + 1)^ + ^:-l+2a-;> 



or 



^1 



Remark 9.1 Here we observe again the cancelation of the leading exponential terms and stress again 
the importance of subdominant exponentials for the calculation of the asymptotics of eigenfunctions 
on all intervals {qj,qj+i). 



■ (j) , I 

We finish this section by computing the coefficients D)^ from (52). Under assumptions (67), we 



have 



D 



(1) 



[1 + 



(1) 



D 



(2) 



b-l 



b-a\ 



D 



D 



(3) 



(3) 



-l + £ 



-eh 



b—a 



[l + £ 



b-a\ 



(68) 



D 



(4) 



b—a 



D 



j^{i + £ 

b — a 

-l + £^-^. 



b-a\ 



Remark 9.2 j^There exist trigonometric polynomials / satisfying assumptions of figure 16, i.e. having 
two local minima and two local maxima satisfying inequalities (64), or (65), or (66), or (67). Indeed, 
one should take any Morse C°° function /o with two local minima and two local maxima satisfying 
the inequalities, say, 

/o(gi) < /ofe) < k{qi) < /ofe); 2[/o(g4) - /ofe)] < 3[/o(g3) - /ofe)]; 2[/o(g4) - Mqi)] < /o(92) - /ofe), 

(69) 

that are, up to shift and rescaling, correspond to the conjunction of (67) and (63). Then the Fourier 
series of /o will converge to /o uniformly together with all derivatives, and an n-th partial sum /„ of 
that Fourier series for sufficiently large n will have critical points and critical values arbitrarily close 
to those of /. Since our conditions (69) are open, will satisfy them for n large enough. With a 
little more work one can produce a trigonometric polynomial with exactly prescribed critical points 
and critical values. Alternatively, one can generate examples of trigonometric polynomials satisfying 
(67) and (63) using a computer algebra system. 



^The material contained in this remark was explained to the author by Prof. A.Gabrielov. 
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A Useful formulas 



For ^ > and E > and E — )• 0+ we have the fohowing asymptotics of various integrals: 

arccosh (A/Ve) = Ln2A- ^Ln E- ^ + o{E) (70) 



cosh'^ tdt 



arccosh (^s=) 



-arccosh {A/yE) 

-^ + \ + TQ^ - sarccosh {A/^/E) + o{E) 



-A^E-i 



4E 32 

k 



farccosh (A/VE) + o{E^) 



2(k-2) 



8(fc-4) 



+ oiE^ 



iik = Q 
iik = 2 
if /c = 4 

if /c = 1, 3 or > 5 
(71) 



f -4 + i + + ^arccosh (^//E) + o{E) if k = 

/ sinh^ t cosh^ tdt=i + M + |arccosh (A/VE) + o(^0) if A; = 2 

[ i?--t (-4?? + # + + if = 1 or > 3 



The following formulae are simple integration by parts used in sec|4j 

du 1 u f udu 1 



(^2 _ ^)3/2 E ./^J^ ' J {U^- ^)3/2 V^J^ 



du I u 1 u'^ f udu 



{y2-E)l E^ Vn2 - E E^3{u'^-E)l ' J {u"^ - E)l 3(n2-£;)i 



2 



(72) 



(73) 



(A; - 1) / dn, A; > 2 (74) 



(75) 



u'^du u^ ^ k — I , ^ , , ,„ , 

5 = 3 + / wdu, k>2 76 

(^2-^)5 3(tx2-S)i 3 J (u^-E)l 
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